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MATHEMATICAL  MODELLING  OF  THE  RAE  LIGHTWEIGHT  FLEXIBLE  SOLAR  ARRAY 

by 

R.  N.  A.  Plimmer 


SUMMARY 

The  paper  describes  the  development  of  a mathematical  model  for  the 
attitude  dynamics  of  a spacecraft  equipped  with  a lightweight  flexible  solar 
array  of  RAE  design.  The  theory  has  been  developed  using  a continuous  mechanics 
approach  and  a computer  programme  prepared  to  generate  the  lateral  bending  modes 
of  a spacecraft  comprising  a rigid  central  structure  carrying  a pair  of  solar 
arrays  symmetrically  situated  about  the  central  body.  Furthermore  the  programme 
will  generate  the  effective  inertia  and  mass  as  a function  of  the  forcing 
frequency.  These  are  then  formulated  in  the  form  of  transfer  functions  which 
are  more  convenient  for  control  problem  analysis.  

Two  models  are  considered.  The  first  model  assumes  that  the  array  cross 
members  are  rigid  whilst  the  second  takes  account  of  the  flexibility  of  the 
cross  members. 


Paper  presented  at  the  VI  IFAC  Symposium  "Control  in  Space",  Teakhkadzor , 
Armenia , USSR,  26-31  August  1974. 
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I INI ADDUCTION 


As  the  power  requirements  for  future  generations  of  communication  and  scientific 

satellites  continue  to  escalate,  radically  different  design  departures  must  be  made 

from  the  early  spinning  drvm  type  of  satellite  with  body  mounted  solar  arrays. 

These  designs  have  already  included  rigid  fold-up  arrays  (Skylab) , semi-rigid  fold- 

19  arrays  (Boeing  Mars  mission),  flexible  roll-up  arrays  (FRUSA) , flexible  fold-up 

arrays  using  BISTEMS  (CTS)  and  flexible  fold-up  arrays  using  a central  telescopic 

mast  (X4)  and  are  usually  sun  orientated  from  a three  axis  stabilised  spacecraft. 

Designs  embodying  flexible  panels  have  been  shown  [l]  to  be-  superior  to  rigid  and 

semi-rigid  types  in  terms  of  power-to -weight  ratio,  stowed  volime,  adaptability  and 

development  potential.  The  estimated  end-of-life  power  to  weight  ratios  of  a two 

paddle  2kW  array  of  the  thinne:;*  available  soi.r  cells,  for  instance,  complete  with 

orientation  and  power  transfer  mechanisms,  range  from  20  W kg"1  for  a rigid  fold-up 

.1 

type  to  43  W kg  for  a flexible  fold-up  type  deployed  by  a pneumatic  telescopic 
aast.  At  the  RAE  (Space  Department)  effort  has  been  concentrated  into  the  develop- 
ment of  the  latter  type  of  lightweight  flexible  array  [2,3]. 

Fig.l  is  an  artist's  impression  of  the  RAE  type  array  fully  oeployed,  as 
originally  proposed  for  a spacecraft  capable  of  demonstrating  orbital  manoeuvres 
by  ’metric  propulsion  whilst  Fig. 2 illustrates  the  main  features  of  the  RAE  proto- 
type wi,i.ch  has  been  developed  over  the  past  sb  /ears  and  has  now  satisfactorily 
completed  qualification  tests  for  geostationary  orbit. 


Ikie  to  the  lightweight  structure  of  the  array,  the  design  of  an  attitude  control 
system  for  a spacecraft  of  the  type  illustrated  in  Figs.l  and  2 has  to  take  into 
account  the  flexible  nature  of  the  array  since  the  natural  vibration  modes  of  the 
spacecraft  may  be  excited  by  the  control  torques  and  forces  applied  to  the  space- 
craft for  the  purpose  of  attitude  control,  station  keeping  and  orbit  manoeuvres. 

The  present  paper  describes  a method  which  has  been  employed  to  assess  the  various 
modes  of  vibration.. 


2 BRIEF  DESCRIPTION  OF  THE  RAE  ARRAY 

Each  of  the  two  arms  of  the  deployed  array  consist  of  a number  of  panels  (six 
pairs  in  the  current  design)  of  thin  (100  pm)  silicon  cells  mounted  on  flexible 
Kaptan  sheets  supported  by -aluminium  honeycomb  cross  members  attached  to  a pneumat- 
ically deployed  telescopic  mast.  The  mast  is  made  up  of  a number  of  thin  walled 
aluninium  alloy  tubes  of  progressively  decreasing  diameter  and  wall  thickness, 
adjacent  sections  of  tube  being  locked  together  at  the  overlaps  to  give  a fairly 
rigid  connection.  The  cross  memberr  are  attached  to  the  outboard  ends  of  the 
tubular  sections  and  the  panels  are  held  under  tension  by  springs  connected  to 
adjacent  cross  members.  During  launch  each  arm  of  the  array  is  folded  concertina 
fashion  and  stowed  in  an  aluninium  honeycomb  compartment,  the  cross  member  on  the 
outermost  tube  forming  the  cover  of  the  compartment.  Once  deployed  in  orbit,  the 
present  design  of  array  cannot  be  retracted. 
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3 HTF.1J-ING  ASStMTIONS  FOR  IMS  ARRAY 

Although  such  arrays  are  frequently  Investigated  nowadays  by  using  the  Method 
of  finite  elements  in  the  final  phase  of  analysis,  this  can  he  expensive  and  time 
cnn«»ing  in  the  initial  phase  of  design.  The  approach  enployed  in  this  paper  has 
been  to  use  the  traditional  methods  of  continuum  mechanics  rather  than  a discrete 
method  since  the  geometry  of  the  system  is  relatively  simple  and  the  method  can 
give  a quick  insight  into  the  problems  involved.  This  method  has  also  been 
successfully  used  by  ftighes  [4,5,6]. 

Although  torsional  and  in-plane  vibrations  will  occur,  initial  estimates  shewed 
that  these  modes  of  vibration  are  at  much  higher  frequencies  than  the  dominating 
out-of-plane  modes.  Consequently  only  the  lateral  flapping  modes  of  vibration  are 
considered  in  this  paper. 

In  order  to  model  the  system,  certain  idealisations  have  been  made:- 

(i)  the  sections  of  the  mast  have  been  treated  as  uniform  cylindrical  beams  jf 
specified  stiffness  and  density  under  compression  due  to  the  panels  and  the 
overlap  portions  of  the  mast  neglected  apart  from  their  mass  contribution. 
Rotary  inertia  and  shear  are  neglected. 

(ii)  the  panels  have  been  idealised  as  uniform  membranes  of  specified  density 
under  uniform  tension  but  with  negligible  flexural  rigidity. 

(iii)  the  cross  members  have  been  assumed  rigid  in  the  first  instance  so  that  they 
can  be  considered  as  point  masses  possessing  inertia.  Flexible  cross  members 
are  considered  later  in  the  paper. 

(iv}  the  central  body  is  assuned  to  be  rigid  with  the  arrays  being  deployed 
symmetrically  about  the  centre  of  mass. 

(v)  small  linear  displacements  and  zero  energy  dissipation  is  assumed,  i.e.  zero 
clasping. 

A ANALYSIS  OF  THE  LATERAL  MOTION 

Fig. 3 shows  the  coordinate  system  used  to  describe  the  mode  shape  of  the 
idealised  array  at  the  ith  section  (i  * 1 to  N) . [The  index  i will  be  omitted 
*4»erever  possible  in  the  text  to  avoid  unnecessary  symbolic  complexities.]  The 
quantity  u(y)  denotes  the  mast  displacement  from  a fixed  reference  plane  whilst 
v(y)  represents  the  corresponding  panel  displacement. 

Since  the  panels  axe  under  tension  T and  each  mast  section  under  the  same 
compression  force,  the  equations  governing  their  motion  are  respectively 

Tv"  ♦ pv  • 0 (1) 

and 


ou""  ♦ Tu"  ♦ pu  * 0 


(2) 


(using  the  notation  prime  to  signify  partial  differentiation  with  respect  to  y 
and  a d-'t  to  denote  partial  differentiation  with  respect  to  tire) 

where  p * miss  per  unit  length  of  the  panel 
p - mass  per  unit  length  of  the  mast 
E * Young’s  modulus  of  the  mast 
I * area  moment  of  inertia  of  mast  » 
a - El 

r j » outer  radius  of  the  mast  » r^  + h 
r 2 * inner  radius  of  the  mast 
h * wall  thickness  of  the  mast 
t * length  of  the  mast  section. 

Thus,  when  the  system  is  vibrating  with  sinusoidal  frequency  p , the  solution  to 
the  panel  equation  (1)  is 


vfv)  - u(0  sin  W ♦ »(0)  sin  k (t  - y) 
' s in  Tt 


(3) 


where  u(0)  and  u(l)  denote  the  mast  displacement  at  the  beginning  and  end  of 
the  section  respectively  and  k • p JljT . Similarly,  it  may  be  shown  [7]  that  the 
mast  displacement  is  given  by 

(«2  ♦ «2)u(y)  - («2  cos  ey  * t2  cosh  «y)uQ  ♦ ^ 

. I (cosh  v - «<»  .,»b  • i (aiair  - ^)v0  („ 

where  M • au"  is  the  mast  moment  and  V » -M’  . The  mast  shear  S is  given  by 
S * V - Tu’  , i.e.  it  is  dependent  upon  the  compression  force  and  the  quantities 
fi,  7,  4 and  « are  defined  by 


(S) 


Consequently,  using  equation  (4),  we  may  derive  the  transition  matrix  Un  - [u  ] 


relating  the  quantities  u,  u’,  M 
namely 


and  V at  the  two  ends  of  the  mast  section, 


u 

u 

u’ 

- \f* 

u’ 

M 

M 

V 

l 

_V_ 

where  the  elements  of  (e2  ♦ a2)^  are  given  by 
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»2c h ♦ A2c 

ilch-c) 

1 / s sh  \ 

a v ' 4 J 

it  (<sh  - As) 

1 2ch  ♦ A2c  ~ (Ash  ♦ Cfl) 

s <c  - ch> 

oS2«2(ch  - c) 

nAt(fsh  - As)  A2ch  * <2c 

— «s  - *r,h 

-aA2«2(Ash  ♦ «s) 

-<»A2<2(ch  - c)  -A“sh  ♦ e^s 

2 2 

A*ch  + * c 

the  notation  s ■ 

sin  «t,  c ■ cos  <1,  sh  * sinh  At  and 

ch  ■ cosh  St  . 

Consider,  now,  the  forces  and  moments  acting  at  a cross  member.  These  are  shown 
in  Fig. 4 at  the  junction  between  the  ith  and  (i*l)th  sections  where  the  suffix 
notation  ♦ and  - has  been  used  to  denote  'just  before'  and  'just  after’  the 
junction.  The  total  lateral  forces  F+  are  equal  to  the  sum  of  the  corresponding 
shears  in  the  mast  and  the  components  Tvj  arising  fiom  the  tension  in  the 
panels.  Thus 


Ft  - s±  ♦ TV;  - v±  - Tu;  ♦ Tv^  . (7) 

The  moments  acting  on  the  ith  cross  member  arise  purely  from  the  neighbouring 
moments  M+  K M^j(O)  and  M_  ■ M^(l)  existing  in  the  mast  sections.  Thus  if  wo 
now  introduce  the  state  variable  X ■ (u,u',M,F)  we  may  write 


— - 

— “ 

- - 

— - 

u 

U 

u 

u 

u' 

- U+ 

u' 

and 

u' 

- if 

u’ 

M 

M 

M 

M 

_V_ 

0 

_F_ 

0 

_F  _ 

t 

-v_ 

(8,9) 


where  the  matrices  IT*  are  derived  by  using  equations  (3),  (6)  and  (7)  to  be 


and 


1 

0 

0 


T'(cos  kt  - u^) 


0 

1 

0 

T - T'u 


12 


0 

0 

1 

T'u 


13 


(10) 


IT 


1 

0 

0 

T’(cos  kl  - Ujj) 


0 

1 

0 


'u  - T 
12 


0 

0 

0 

4 


(11) 


where  T'  • Tk/sin  kt  and  A • 1 + T'u^  . Furthermore,  the  equations  of  motion 

for  the  cross  member  of  mass  and  moment  of  inertia  mc  and  I are 

c 
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F+  - F_  - * mru 
and 

M+  - M.  - Icu’ 

giving  the  transition  nvatrix  if  defined  by 
X+  - lfx_ 

- r 1 0 0 

0 1 0 

0 -ICP2  1 

-mcp^  0 0 


(12) 


Hence,  using  the  transition  matrices  defined  by  equations  (6),  (8),  (9)  and  (12) 
the  overall  transition  matrix  A for  a single  array  relating  the  state  at 

the  end  of  the  mast  to  that  corresponding  to  the  connection  of  the  mast  to  the 
satellite,  i.e. 

“ "o-  i (13) 

is 

N 

* ■ TT  (“i «; "?  "0  • iM> 

i-i 

i 

We  will  denote  the  elements  of  A by  a^  (i,j  • 1 to  4). 

5 MODES  OF  VIBRATION 

In  view  of  the  symmetry  assumptions  made  concerning  the  deployed  state  of  the 
two  arrays,  the  lateral  deflections  can  be  considered  as  the  superposition  of  two 
classes  of  normal  mode  responses.  The  first  class  consists  of  those  modes  in  which 
the  two  masts  bend  in  opposite  directions  (i.e.  antisyrmietrical  deflections)  so 
that  there  is  rotation  of  the  main  body  but  no  translation.  The  second  class  is 
conprised  of  the  modes  in  which  the  tvo  masts  bend  in  the  same  direction  (i.e. 
symmetrical  deflections)  and  there  is  thai  translation  of  the  central  body  but  no 
rotation.  These  will  be  referred  to  as  the  rotational  and  translational  modes 
respectively.  For  the  attitude  control  system  it  is  only  the  rotational  modes 
which  are  of  interest  since  it  is  these  which  affect  the  pointing  accuracy  of  the 
satellite.  Another  class  of  modes  will  be  called  the  cantilever  modes  and  refers 
to  the  modes  of  vibration  occuring  when  the  satellite  is  rigidly  fixed,  i.e.  of 
infinite  inertia  and  mass. 

In  order  to  determine  the  various  modes  of  vibration  the  appropriate  boundary 
conditions  nust  be  applied  to  the  transition  equations  (13).  These  will  now  be 
considered  in  turr . 
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6 ROTATIONAL  AND  rAOTI LEVER  HOPES 

If  the  rotation  of  the  satellite  main  body,  of  mas.'  and  inertia  1^  , is 
denoted  by  8 , the  arm  length  from  the  centre  of  mass  to  the  array  attachment 
point  by  d and  the  applied  couple  about  the  centre  of  mass  by  C then  the 
boundary  conditions  at  0_  are 

u - <#  , u'  - 8 , It>  - 2M  ♦ 2Fd  ♦ C 

and  at  N+  , M » F ■ 0 . From  these  conditions  applied  to  equation  (13)  it 
follows,  after  some  straightforward  manipulations,  that  the  effective  inertia 


Ie  - C/8 

of  the  satellite  is 

ZNq 

VP>  ' Jb 

R 

(IS) 

where 

NR  " 

^31  * a3Z 

a34  ■ ^33 

(16) 

^41  ♦ a42 

a44  ‘ ^43 

and 

°c  ■ 

a33  8 34 

a43  a44 

(17) 

are  independent  of  the  central  body  inertia. 

Now  the  natural  rotational  modes  of  vibration  are  those  obtained  when  the 
applied  couple  is  zero,  i.e.  ■ 0 . Also,  the  cantilever  frequencies  Hn 

occur  when  the  satellite  rotation  is  zero  under  an  applied  couple,  i.e.  is 

infinite  or  Dc^n)  ■ 0 . Thus  by  plotting  Ie  j and  Dc  as  a function  cf  p the 
rotational  and  cantilever  frequencies  may  easily  be  determined.  Furthermore, 
equation  (15)  shows  that  the  dependence  of  the  rotational  frequencies  on  the 
central  body  inertia  may  be  deduced  simply  by  a shift  in  the  Ig  ordinate  axis. 

It  ma>  also  be  shown  that  in  the  limit  as  p ♦ 0 , Ie(p)  tends  to  the  overall 
inertia  Ir  of  the  satellite  considered  as  a rigid  body. 

7 TRANSLATIONAL  MODES  OF  VTBRATICN 

Denoting  the  lateral  displacement  of  the  satellite  centre  of  mass  by  Z and 
the  applied  force  by  P then  the  boundary  conditions  at  0 are 

u - Z , u'  - 0 , M^Z  » 2F  ♦ P 

whilst  at  N+  , M ■ F ■ 0 . Again,  applying  these  conditions  to  equation  (13) 
leads  to  an  effective  overall  mass  Mg  • P/Z  given  by 


mmmwm  mmm******-** 


.yt*m  w-  3*wn»««wssR^  .- 


where 


a33  a31 

a43  "l41 


(19) 


which,  as  to  be  expected,  is  independent  of  the  arm  distance  d . Furthermore,  in 
the  limit  as  the  vibrational  frequency  tends  to  zero,  Mp (p)  tends  to  the  overall 
rigid  mass  Mr  of  the  satellite.  The  natural  translational  frequencies  of  vibra- 
tion u/  are  those  existing  when  the  applied  force  is  zero,  i.e.  - 0 end, 

of  course,  M <P  ) is  infinite  at  the  cantilever  frequencies.  As  with  the 
rotational  natural  frequencies,  the  translational  natural  frequencies  niay  be 
determined  by  plotting  Mg  as  a function  of  p and  noting  the  zeros.  Further- 
more, as  can  be  seen  from  equation  (18),  the  effect  of  the  central  body  mass  on 
the  frequencies  may  be  found  again  by  a corresponding  shift  in  the  ordinate  axis 


8 TRANSFER  FUNCTION  RFJ’RESENTATION  OF  THE  SATELLITE 


Knowing  the  natural  modes  of  vibration  we  may  thus  express  the  response  to  the 
applied  forces  and  couples  in  the  form  of  transfer  functions,  namely 


Z 

P 


MrS 


TT 

r . , *2  i 

°J(1  + ifn> 

TT 

II 

r . , s*  i 

(20) 


and 


TTfi  ♦ -y-*! 1 

1 j>L 


O*  — — r r 2 

r JT  1 ♦ — 

• 'L  “n(1  + ilrR> 


(21) 


where  s • ip  aru.  structural  damping  has  been  introduced  in  the  customary  manner 
[8J,  i.e.  the  rigid  inertia  and  mass  have  modification  factors  due  to  the 
flexibility  of  the  array  so  that  the  control  block  diagram  is  relatively  simple 
for  the  particular  configuration  chosen  here  (although  it  may  be  readily 
generalised)  and  only  requires  knowledge  of  the  natural  frequencies.  In  practice 
anly  sufficient  factors  in  equations  (20)  and  (21)  would  be  retained  to  cover  the 
bandwidth  of  the  control  system. 


This  form  is  in  contrast  to  the  form  employing  nodal  gains  [4,5,6]  which  require 
the  actual  mode  shapes  to  be  calculated  and  use  of  thei-  orthogonality  properties. 
The  orthogonality  conditions  for  the  present  model  can  easily  be  deduced  by  consid- 
eration of  the  total  kinitic  energy  cf  the  system  and  using  tiie  'conjugate 
property'  [9].  In  essence,  the  transfer  functions  depends  upon  a mixed  representa- 
tion of  the  constrained  and  unconstrained  modes  of  [4,6]. 
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9 COMPUTER  PRPgjWE 

In  order  to  obtain  the  effective  inertia  and  mass  as  a function  of  frequency  a 
computer  programme  was  written  to  perform  the  matrix  nultiplications  indicated  by 
equations  (14),  (IS)  and  (18)  as  a function  of  the  basic  physical  parameters  (T, 

|i,  E,  p,  Tj,  h,  t)^,  M^,  1^  and  N . The  program*  was  also  organised  to  iterate 

to  the  zeros  of  Ie,  Mg  and  Dc  in  order  to  obtain  the  rotational,  translational 

and  cantilever  frequencies  and,  if  necessary,  to  print  out  the  corresponding  mode 
shapes.  (Normalised  to  have  unity  displacement  at  the  tip  of  the  mast.) 

As  an  illustration  of  some  of  the  results  obtained  Figs.S  and  6 show  the  effect- 
ive inertia  and  mass  obtained  for  tne  array  having  the  physical  characteristics 
shown  in  Table  1.  It  will  be  seen  that  the  rigid  inertia  of  the  two  arrays  is 
390  kg  m2  whilst  the  lowest  rotational  frequency  for  a central  body  inertia  of,  say, 
300  kg  m 2 is  0.44  Hz  and  decreasing  to  the  cantilever  frequency  of  0.3  Hz  as  the 

central  body  inertia  is  increased.  The  behaviour  at  0.9  Hz  and  1.8  Hz  is  not  shown 

in  detail  since  there  are  other  very  close  rotational  frequencies  corresponding  to 
panel  excitation.  These  rotational  frequencies,  and  their  variation  with  mast 
stiffness,  are  shown  in  Fig. 7 for  1^  ■ 300  kg  m2  . 

I 

Physical  Characteristics  of  Single  Array  (Mass  - 9.20  kg) 


Section 

t 

(m) 

F 

(kg/m) 

Ft 

(kg) 

rl 

(nw) 

h 

(mm) 

9 

(kg/m) 

fit 

(kg) 

rac 

(kg; 

1 

1.295 

0.9594 

1.243 

25.40 

0.4572 

0.2180 

0.2825 

0.0P10 

2 

1.29S 

0.9594 

1.243 

25.40 

0.4572 

0.2180 

0.2825 

0.0288 

3 

1.295 

0.9594 

1.243 

22.23 

0.4572 

0.1766 

0.2288 

0.0288 

4 

1.295 

0.9594 

1.243 

20.64 

0.4572 

0.1623 

0.2102 

0.0288 

5 

1.295 

0.9594 

1.243 

19.05 

0.3310 

0.1250 

0.1619 

0.0288 

6 

1.295 

0.9594 

1.243 

17.46 

0.3810 

0.1143 

0.1481 

0.2268 

Total 

7.770 

7.458 

1.314 

0.424 

Fanel  tension  5N,  Young's  modulus  of  mast  7.03  * 1010  N in'2,  width  1.37  m 


10  EXTENSION  OF  THE  fDDEL  FOR  FLEXING  CROSS  MEMBERS 

For  weigh'  reduction  purposes,  the  possibility  of  replacing  the  stiff  cross 
members  by  relatively  flexible  steel  hyperdermic  tubing  is  being  seriously 
considered.  In  order  to  determine  the  consequent  change  in  flexibility  of  the 
array,  another  model  has  been  developed  which  will  now  be  very  briefly  described. 
The  coordinate  system  is  similar  to  that  of  Fig. 3 except  that  the  panel  motion 
v(x,y)  is  now  a function  of  x and  y whilst  w(x)  represents  the  displacement 
of  the  cross  member  from  its  undeflccted  position  and  is  represented  by  a super- 
position of  clamped- free  beam  modes  ^(x)  [lO],  i.e.  w(x)  » Eq^.(x)  with 

*i(0)  * V- (0)  - 0 . The  igj  satisfy,  in  particular,  the  equations  [ll] 
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h *i  " C<i3  • l*i  ' f~  * ^ “ *}«y 

^ C f. 


(22) 


where  the  suffice  c on  the  integral  denotes  integration  over  the  length  of  the 
cross  mother  and  the  quantities  «.  and  fi^/2  are  tabulated  in  [lOj. 

In  order  to  develop  the  equations  of  motion  recourse  may  be  made  to  Hamilton's 
principle,  namely 

«/(!.♦  WH’t  - 0 (23) 

In  which  4 denotes  a virtual  variation,  l is  the  array  Lagrangion  and  w is 
the  mrk  done  on  the  array,  e.g.  by  the  forces  and  r*nents  at  the  root  of  the  array. 
Here  the  Lagrangian  may  be  written  as 

1 " / p"?dy  * } MV^dxdy  *f£  („  ♦ 

L-  p c 

u,  2dy  ♦ J Tu'^dy  - /l  V'2dxdy  - / «c  (S^f  ^ 

* P c W/ 

where  £ indicates  station  over  all  sections  of  the  array  and  ac  is  the 
bending  stiffness  of  the  cross  member. 

Carrying  eut  the  variation  in  equation  (21)  and  (24)  lead,  again  to  the  partial 
differential  equations  (1)  and  (2)  with  solutions  (4)  and 


7 


- / au' 


(24) 


. ,,w  sgvg . v0w 


(2S) 


V,her0  VQ(x)  - u(0)  ♦ ^ ^(x) 

Vt(x)  - u(lj  ♦ Xqjy^x} 

together  with  the  boundary  conditions  at  a cross  member 

S+  ‘ S-  + /(«  V'ldj<  ’ /(e  v’).dx  + P2"»C(u  * la’q^  . o 


and 


/(!  ^.dx  - J(l  v'j^.dx  - Kiq.  ♦ pV(uo-  ♦ q.)  - 0 


(26) 
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(i  “ 1,2  ...) 


K,  - «CdJc  . Using  equation  (25),  we  may  evaluate  the  quantities 


r* 


J / (i v')+dj<  " 4 f28) 

"i  ■ /(?  v’)dx  ■ T'  [(ut  * cos  kl  “ (“0  + ^ (29) 

J;  C 

l /(I  v’)/idJC  ■ r [(vi * ql)  * (vi  * q?) cos  kt]  a <?  . <“> 

«£ 

I / (f  v']/idx  " T’  [(vi  * qi)  cos  kt  ■ (vi  * qi)]  E Gi  (31) 

f and  equations  (26)  and  (27)  become  respectively 

j F*  - F*  ♦ p2mC (u  ♦ Sq.ap  - 0 (32) 

► 

f and 

j Gi  * Gi  ■ Kiqi  * P^fua!  ♦ - 0 (33) 


(i  - 1.2,  ..0 

where  F « S ♦ . Furthermore,  we  have  the  equations  of  continuity  at  a cross 

member,  namely 

u*  - u'  , u'*  - u'“  , M*  - M*  , qi  " qi  (34) 

where  we  have  assumed  that  the  moments  of  inertia  of  the  cross  members  are 
negligible. 

If  we  terminate  the  expansion  for  w(x)  by  retaining  only  n terms  then 
equations  (32),  (33)  and  (34)  enable  us  to  obtain  the  transition  matrix  ll*"  at  a 
cross  member  relating  the  state  variable  Y ■ (q^  ...  q^  u u'  M F Cj  ...  G^)  of 
order  2n  ♦ 4 at  ' ‘ust  before'  to  'just  after'  the  junction,  i.e. 


I y.  » ify 

! 

j where  the  ele.uents  of  Uc  are  given  by 


/ 
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and  1 . 0 denote  unity  and  null  matrices  of  order  n respectively, 
n’  n 

Using  the  relation  S - V - Tu'  , equation  (6)  may  be  modified  to  be 


- — 

-*  — 

u 

U 

u* 

. u*m 

u' 

M 

M 

_S 

S 

t 

L- 

(35) 


where 


.*m 


11° 


0 

0 

0 

irIta2i 


Tu 

Tu 

Tu 


Tu44  - T(u 
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l24 

34 


22  + Tu24^ 


with  elements  Vy  . 

Consider  now  equations  (28),  (29),  (30),  (31),  (35)  and  the  relation 
pQ  t • (S  ♦ ( . These  give  in  ♦ 8 equations  relating  the  2(2n  ♦ 6) 

quantities  (q^  ...  u u*  MSF^jGj  ...  Gn)Q  t . By  eliminating  the  four 
quantities  (S,  C^)Q  t we  are  left  with  (2n  ♦ 4)  equations  relating  the 
2(2n  ♦ 4)  elements  of  i * i,e‘  a trailsttion  equation  in  the  form 

Vt  ‘ Vo  M 

for  the  ith  section  of  the  mast.  Besides  the  2n  equations  (30)  and  (31),  the 
other  four  equations  are 

Iv14T'aiqi  + (1  * v14T'iut  “ rv14T’  cos  U °iqi  + (vH  + vi4T'  cos  ki)u0 

+ V1 2 “6  + v13^  * v14F0 

rv24T'aiqi  * v24T'ut  + ui  " Zv24T'  cos  **  aiqi  + (v21  * V24T'  cos  ktJu0 

+ v22uO  * V23^0  + v24F0 

JV34T’tfiqi  + v34T’ul  + Ml  ' Xv34T'  ccs  kjl  °iqi  + (v31  + V34T’  ODS  kt)uo 

+ V32U0  + v33*b  + V34F0 

2(v<4  - cos  ki)T'ajq*  ♦ (v44  - cos  kl)T'ut  + Fj 

- r(v44  cos  kt  - l)T'a!q°  ♦ (v4I  ♦ v44T'  cos  kt  - T')^  ♦ v42u£  ♦ v ^ + v44FQ 
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so  that  the  elements  of  TQ  and  Tj  are  determined  explicitly.  Consequently,  the 
overall  transition  matrix  B for  a single  array  relating  the  state  at  the  end  of 
the  mast  to  that  at  the  root,  i.e. 


B - 

Knowing  B , the  effective  mass  and  inertia  for  the  translational  and  rota- 
tional motion  can  be  found  in  exactly  the  same  manner  as  the  earlier  discussion  by 
applying  the  appropriate  boundary  conditions.  In  particular  it  is  to  be  noted  that 
■ 0(i  ■ 1 to  n)  at  0_  and  N+  . A computer  programme  lias  also  been  written 
for  the  flexible  cross  member  model.  Numerical  results  will  be  given  at  the 
symposium. 
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